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Abstract. Essential properties of halo nuclei can be described in terms of a few low-energy constants. 
For neutron halo nuclei, analytical results can be found for wave functions and electromagnetic transition 
matrix-elements in simple but well-adapted models. These wave functions can be used to study nuclear 
reactions; an especially simple and instructive example is Coulomb excitation. A systematic expansion 
in terms of small parameters can be given. We present scaling laws for excitation amplitudes and cross 
sections. The results can be used to analyze experiments like Be Coulomb excitation. They also serve as 
benchmark tests for more involved reaction theories. 
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1 Introduction 

Exotic nuclei are available as secondary beams at many 
radioactive beam facilities around the world. These un- 
stable nuclei are generally weakly bound with few, if any, 
excited states. A well developed method to study halo nu- 
clei is Coulomb excitation. An instructive example is the 
excitation of the l/2~ bound state in n Bc from the l/2 + 
ground state [TPPP]. 

Halo nuclei are a low-energy phenomenon and can be 
described effectively in terms of a few low-energy param- 
eters. The ratios of core size to the sizes of the halo states 
(to be defined in eq. ([T]) below) serve as small expansion 
parameters. Wave functions and matrix elements can be 
given in terms of these parameters. These wave functions 
can be used in reaction models and simple and realistic 
formulae are obtained. 

In a recent paper the B(E1; 1/2+ — > l/2~) strength 
for 11 Be has been determined from intermediate energy 
Coulomb excitation measurements [1] . In order to analyse 
such kind of data in terms of electromagnetic matrix ele- 
ments it is certainly necessary to use rather sophisticated 
codes which take higher-order electromagnetic and nuclear 
effects into account. These codes can be checked and val- 
idated by comparing their results to limiting cases where 
analytical results can be obtained, e.g., pure Coulomb ex- 
citation that shows some rather simple features. This can 
provide a useful guide for more sophisticated approaches. 
It is the purpose of this paper to provide such analytical 



Send offprint requests to: S.Typel 



results. In the theoretical analysis of [T], quite a compli- 
cated XCDCC method was used, which is required for 
a quantitative description of the experimental data, but 
it may tend to obscure the understanding of the simple 
physical mechanism of Coulomb excitation. 11 Be is an 
archetype of a halo nucleus with a 10 Be core and a sin- 
gle halo neutron in the 2si/2 state. There is a strong El 
transition to the l/2~ bound state, itself a p- wave neutron 
halo state. This dipole transition was previously studied 
by Coulomb excitation at GANIL, RIKEN, and MSU [1 

EE]. 

There are two somewhat separated questions: the influ- 
ence of nuclear excitation in grazing collisions and higher- 
order electromagnetic effects in distant collisions where 
nuclear interaction effects can safely be neglected. We deal 
here with the second question: higher-order electromag- 
netic excitation. This problem was studied for the n Be 
case in the mid-nineties by two groups [MS]- In view of 
recent advances in the description of the electromagnetic 
properties of halo nuclei, see e.g. [71I51IMTTJ] . it seems ap- 
propriate to update this work using the recent analytical 
results for halo wave functions and present scaling laws 
for Coulomb excitation of halo nuclei. Effective field the- 
ory methods are also applied successfully to halo nuclei, 
see [TT1IT5]. 

In section 2 we give the main theoretical formulae for 
our model of the Coulomb excitation of neutron halo nu- 
clei from s- to p-states. In section 3 scaling rules are dis- 
cussed. Then we give an application to the case of 11 Be. 
This can serve as a benchmark for more involved studies 
like XCDCC or time-dependent approaches [13] • Conclu- 
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sions are given in section 4. A preliminary account of part 
of this work was published in the proceedings of the con- 
ference 'Nuclear Physics in Astrophysics III' in Dresden, 
March 2007 M- 



2 Analytically solvable model for Coulomb 
excitation of neutron halo nuclei 

A few low-energy parameters are sufficient to characterize 
halo nuclei. Let us consider the single-particle excitation 
of a neutron from a ground state i = to a bound excited 
state i = I with neutron separation energies Ei > 0. (For 
n Be we have E = 504 keV and E 1 = 184 keV.) The size 
of the single-particle wave functions is determined by the 
bound-state constants qi = ^/2jTEl/h with the reduced 
mass fj, = m n m c /(m n + m c ) of the nucleon+core system. 
With the radius R of the core beyond which the nuclear 
interaction is assumed to vanish, it is possible to introduce 
dimensionless parameters 



qiR 



(1) 



that are a measure for the ratio of the core size to the 
size of the neutron wave functions. These parameters are 
small for halo nuclei, see, e.g., OH], and can be used as 
convenient expansion parameters. In this spirit, we use 
this single-particle model for the jf = 1/2+ and jj = 
l/2~ states in order to evaluate the matrix element for 
the El electromagnetic excitation. It is dominated by the 
exterior contributions. 



2.1 Halo wave functions and matrix elements 

The single-particle wave functions are given by 



*i(r) 



Mr). 



y hs (?) 



(2) 



for the ground state (i = 0) and the excited state (i = 
1). The angular dependence is described by the spinor 
spherical harmonics y 1 -^ (s — 1/2). We use the results 
and notation of [8] assuming zero spin of the core. The 
radial wave functions of the two states are given in the 
exterior region (r > R) by 



/o(r) = C q rh { p {iq Q r) = -C exp(-g r) 



and 



/i(r) = Ciiqirh^ (%r) 



(3) 
(4) 



= -Ci exp(-gir) H 

V Qif 

with Hankel functions hi of imaginary argument. Both 
states are halo states and the normalization cons tants are 
given in the halo limit by Co = \/2qo and C\ = y / 2q 2 R/3, 
respectively [8]. 



The parameters qo and q± of the shallow bound states 
in a halo nucleus are closely related to the scattering length 
ai and the effective range parameter ri for partial waves 
I > in the effective range expansion k 2l+1 cot(<5;) = 
-l/ai+nk 2 /2+. ... The S matrix 5/ = [cot(Si)+i]/[cot(6i)- 
i] has a pole at k = iks, i.e. cot[<5/(fcs)] = i for a bound 
state with ks = qt- This gives the desired relation 



(-l) i+ V +1 = -^-^f 



(5) 



Now there is a difference between I = and I > 0. For 
I = and small go we have in lowest order the relation 
q = 1 /do and the effective-range term is a small cor- 
rection. For / > the left hand side of (0 is smaller 
than the first two terms individually on the right hand 
side and we have an enhanced value (as compared to "di- 
mensional considerations") of the scattering length that 
is given in lowest order by ai = —2/(riqf) ("fine tuned"). 
In [7] the scattering length in the I — 1, j — 1/2 chan- 
nel of 11 Be was determined to be 457(67, —66) fm 3 . With 
qi = 0.0895 fm in this channel we find an effective-range 
parameter r\ — —0.547 fm -1 of "natural order" R^ 1 . 

With the radial wave functions and (fj| we calculate 
the B(E1) value for the l/2 + — > l/2~-transition as well 
as the higher-order effects in electromagnetic excitation. 
We propose this to be a model study and leave the spec- 
troscopic factors equal to one. (They could be adjusted, 
which would result in a quasi-realistic description of the 
n Be system for our purpose.) The B(E1) value is given 

by 



Z {1) e 



R 



(G) 



where — Z c m n /(m n + m c ) is the dipole effective 

charge number and the radial dipole integral is given by 



Roi = I drfl(r)rf (r) 




(7) 



., no do + 271) 

3 (70 + 71) 



In the present paper we are only interested in the halo 
limit, i.e., we keep only the lowest-order term in the ex- 
pansion in 70 and 71. In this approximation we can use 
the exterior radial wave functions, eqs. ([3]) and ([4]) in the 
integral eq. ([7]) down to r = 0. The correction terms from 
using the correct interior wave functions are of higher or- 
der in the expansion parameters 70 and 71. For R — > 
the radial dipole integral goes to zero because the normal- 
ization of the p-wave function tends to zero in this limit. 
Thus R must be kept finite; we choose R — 2.78 fm as 
a realistic value for 11 Be [7J. This value determines the 
asymptotic normalization of the p-wave bound state. We 
find B(E1) = 0.193 e 2 fm 2 , to be compared to the value of 
B(E1) = 0.105(12) e 2 fm 2 obtained from an analysis of the 
GANIL data, see [Tj. This value is consistent with other 
Coulomb dissociation experiments at RIKEN and MSU 
and the value obtained by the Doppler shift attenuation 
method [15] . 



S. Typel, G. Baur: Scaling laws and higher-order effects in Coulomb excitation of neutron halo nuclei 



3 



2.2 Coulomb excitation of neutron halo nuclei and 
scaling laws 

We treat electromagnetic excitation in the semiclassical 
approximation. For high beam energies the classical tra- 
jectory can be taken as a straight line with impact param- 
eter b. In the sudden approximation one can take higher- 
order effects into account in a convenient way. Like the 
related Glauber approximation it is applicable for high 
beam energies and low excitation energies. This is reason- 
ably well fulfilled for 11 Be, even at GANIL energies. In [5 
higher-order effects in the electromagnetic excitation of 
11 Be to the l/2~ bound state were studied using the sud- 
den approximation. In [6] a continuum discretized coupled 
channels approach was adopted. In [16j an analytically 
solvable model for higher-order effects in the electromag- 
netic dissociation of neutron halo nuclei was presented. In 
that work, there was only the transition from an s-wave 
bound state to the continuum. Now we consider the case 
where there is, in addition, a p-wave bound state, as it is 
the case in n Be. In the sudden approximation the dipole 
excitation amplitude from the ground state i to the final 
state / is given by 



^suddc 



(/|exp(— igcoui 



where the transfered momentum is 



Tiq 



Coul — 



vb 



(8) 



(9) 



The impact parameter b is chosen to be in the z-direction; 
this is convenient for the following calculation using polar 
coordinates. The target charge number is denoted by Z 
and v is the beam velocity. The dipole approximation is 

quite well fulfilled, since the dipole effective charge zj^ 
is much larger than the corresponding quadrupole charge. 
The neutron and core masses are denoted by m n and m c 
respectively, the charge of the core is given by Z c . The sud- 
den approximation is applicable for £ = ub/v <C 1 , where 
hui = Eq — Ei is the excitation energy. Even for the com- 
paratively low GANIL energies of about 40 MeV/nucleon 
this is reasonably well fulfilled. The most important inter- 
mediate states are expected to be in the low-energy con- 
tinuum, where the dipole strength has a peak, at around 
1 MeV excitation energy. Thus the adiabaticity condition 
£ -C 1 will also apply to the important intermediate states. 
We note that for the virtual excitation of the high-lying 
giant dipole state the opposite limit £ 3> 1 is realized, this 
would lead to a real polarization potential, see e.g. [T7]. It 
influences the classical trajectory of the projectile and it 
is not important in our context. It is also neglected in our 
model space. The sudden approximation has the advan- 
tage that intermediate states are treated by closure, thus 
one only needs a model for the initial and final states, 
and not for all the intermediate states. In lowest order 
in gcoui the first-order dipole approximation is obtained 
corresponding to a single-photon exchange between target 
and projectile. It is shown in [5] that third-order El ex- 
citation is more important than second-order E1-E2 ex- 
citation. The role of higher multipole excitations in the 



electromagnetic excitation of one-neutron halo nuclei is 
also studied in [18] , 

In the appendix we show that 



^sudden ( 1)^ 

where 



2i 



h{z) = 2 



—h(z) + -I 2 {z) 



arctan(2:) 



and 



h{z) = ^1 + arctan(z) — — 
depending on the dimensionless parameter 

<7Coul 9Coul-R 



(10) 

(11) 
(12) 

(13) 



qo + 9i 7o + 7i 

This quantity is proportional to the usual strength param- 
eter 

(1) _ Ze(f\\M(El)\\i) 



hvb 



(14) 



that characterizes the importance of higher-order effects 
in Coulomb excitation processes. We note that this param- 
eter becomes large in the halo limit (7 — > 0). In contrast 
to x ) the adiabaticity parameter £ = hu>/(vb) with the 
excitation energy hu> = % (q$ — qf)/ (2/i) becomes small in 
the halo limit and the application of the sudden approxi- 
mation {5} is very well justified. 

We can expand the excitation amplitude IjlOp in powers 
of z with the help of the power series of the arctan func- 
tion. This power series expansion converges for \z\ < 1 
and one obtains 



Qsuddcn — 0.1 + CL3 + . . . 

with the lowest-order term 

3i V 3 70+71 
and the next-to-leading order term 



«3 



(— \y i+mi — 



2 1 



70 70 + 471 ^ 3 



15 V 3 70+71 



(15) 



(16) 



(17) 



The even terms vanish due to the parity selection rule. 
Note that in the sudden approximation ||5J) there is no 
change in the magnetic quantum number, since the vector 
<ZCoui is in the ^-direction. 

Following Ref. [5] we can introduce a reduction factor 



r(z) 



| ^sudden | 

kl 2 



1 



2 7o + 47i 2 
Z 

5 70 + 271 



(18) 



In figure[T]the dependence of the reduction factor r on z is 
shown for the present case with 71/70 = \/Ei/E a w 0.604 
by the solid line. We find r(z) w 1 — 0.619 z 2 in low- 
est order in z. It seems of interest to compare the re- 
sults of [S] with the present one. In [S] two different mod- 
els (denoted by / and II) were used. Expanding in the 
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Fig. 1. Reduction factor r{z) = |a su dd cn | 2 / \ a i\ 2 as a function 
of the scaling parameter z as defined in eq. (|13[1 for values of 
(dashed line), 0.604 (solid line) and 1 (dotted line) of the ratio 

7l/70 = y/E 1 /E . 



r 11 (z) = 1 — 12z 2 /5 + . . .. The present higher-order cor- 
rections are smaller than the one found in [S] . In that ref- 
erence rather crude models were used for the first excited 
l/2~-state; those wave functions extended much further 
out than the present wave function which takes the 
angular momentum barrier into account properly. In the 
higher-order terms the importance of the outer region is 
enhanced due to the weighting with a higher power of r. It 
should be recalled that the main aim of [5] was to provide 
an upper limit for higher-order effects in order to under- 
stand the results of [2]. In the meantime, other Coulomb 
excitation experiments were performed Q3I3U1] that clari- 
fied the situation. 

In the general case of a Coulomb excitation from a 
bound s-wave to a bound p-wave neutron halo state, the 
energy E\ of the excited state is limited from below by 
zero (state at the breakup threshold) and from above by 
Eq (zero-energy excitation) . This corresponds to the inter- 
vall [0, 1] for the ratio 71/70- The two limiting cases are 
depicted in figure Q] by dashed and dotted lines. 

The excitation probability is given by 



^(k) — 9 ■ 1 1 ^ 5„ li1nf |a su ddcn| 



2j 

PLo(b) + PnloO) 



(19) 



The lowest-order term 



4 7o (7 + 2 7l ) 2 
27(70 + 7i ) 2 



A-2 

-w (20) 



in the excitation probability is proportional to z 2 . The 
most important higher-order contribution comes from 03. 
Its interference with the lowest-order term a\ leads to the 
next term 



P 



NLO 



(b) 



870(70 + 2 7i)(7o +47i 
135(7o + 71 ) 2 



' 4 = "£ (2D 



in the expansion of P(b) in z, of the order of z 4 
constants A2 and A4 can be written as 



An 



470(70 + 271) 2 / 2??TO r 
27(7o + 7i ) 4 \m c + 



and 



+4 



870(70 + 2 7i)(7o + 4 7i 
135(7o + 71) 6 



m r , 



2rjm, n 



R 2 



The 



(22) 



R 4 (23) 



where the Coulomb parameter is given by 77 = ZZ c e / (hv). 
Total cross sections are obtained by integration over the 
impact parameter, starting from a minimum impact pa- 
rameter 6 m j n . The sudden approximation fails for large im- 
pact parameters, and an adiabatic cut-off 6 max = Tbeam^/w 
with the Lorentz- factor by 7bcam = 1/ v/l ~~ (v/c) 2 has to 
be introduced for the lowest-order result. For the higher- 
order terms this is not necessary, the convergence in b is 
fast enough. We get 



ctlo = 27T 



P LO (b)bdb = 2TrA 2 In ■ 



strength parameter z we get r (z) = 1 — 6z 2 /5 + . . ., and anc j 



cnlo = 2-7T 



P N Lo{b)bdb 



nAi 

mm 



(24) 



(25) 



(A somewhat different method was used in eqs. (16) and 
(17) of [5], where the ^-dependence of the first-order am- 
plitude cii is taken into account. Then there is no need for 
the adiabatic cutoff b max .) 

We note that the strength parameter z is proportional 
to 1/v, i.e. the leading-order term decreases like l/E, the 
next-to- leading-order term like 1/P 2 , where E is the beam 
energy. In the halo limit (7 = 70, 71 — > 0) the probabilities 
scale as Plo oc I/7 and Pnlo oc I/7 3 [z scales as I/7) and 
the excitation probability tends to infinity. However, the 
NLO contribution tends to infinity even faster and higher- 
order corrections become more important. Still, for realis- 
tic values of 7 the higher order effects are quite small (see 
below). It seems intuitively understandable that higher- 
order effects tend to increase with decreasing 7. In fact, 
this increase with decreasing 7 is faster than the increase 
in lowest order, as our analysis shows. 

For the total cross section there is the additional well 
known enhancement In 7 in LO, which is absent in NLO, 
due to the fast convergence in b. 

In [16] and [19] the scaling properties of the Coulomb 
dissociation of halo nuclei were investigated in a related 
approach. The strength parameter defined there corre- 
sponds to the one defined now. Using the present nota- 
tion we found in that case that Plo oc I/7 3 and Pnlo °c 
I/7 5 . The difference to the present case arises because the 
bound-continuum and bound-bound state dipole matrix 
elements show different scaling properties with 7^. 

3 Application to the 11 Be Coulomb excitation 
experiments 

In the case of 11 Be we have neutron separation energies 
of E Q = 504 keV and E x = 184 keV for the ground 
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Table 1. Projectile energy per nucleon -E/A, velocity v, max- 
imum impact parameter b max , maximum parameter z max and 
cross section reduction factor R(v) for the n Be experiments. 





E/A 


V 


^max 




R{v) 




[MeV] 


[c] 


[fm] 






GANIL [TJ 


38.6 


0.2861 


184.1 


0.6549 


0.8902 


GANIL [2] 


43.0 


0.3009 


194.6 


0.6227 


0.9025 


MSU g] 


60.0 


0.3507 


230.9 


0.5343 


0.9321 


RIKEN [3] 


64.0 


0.3611 


238.8 


0.5189 


0.9366 



and excited bound state. This leads to an excitation en- 
ergy of hui = 320 keV for the l/2~ state. With a radius 
R = 2.78 fm one obtains the numerical values 70 = 0.4116 
and 71 = 0.2487, respectively, for the dimensional scaling 
parameters that appear in the excitation amplitude. De- 
pending on the projectile energy, the projectile velocities 
v and the maximum impact parameters 6 max = Jbc&mV / w 
are found for the conditions of the experiments at GANIL 
[TP]. RIKEN [3J and MSU g]. They are given in ta- 
ble [JJ These adiabatic cutoff radii are impressively large. 
It means that Coulomb excitation extends really far out 
and there is an amply large region where the interac- 
tion is purely electromagnetic, and our approach works 
very transparently. Of course, there is a Coulomb and nu- 
clear interference zone, quite moderate in extension, close 
to the minimum impact parameter 6 m i n = 1.2(11 1//3 + 
208 1 / 3 ) fm = 9.78 fm which is not considered in our ap- 
proach. 

The maximum value of the dimensionless strength pa- 



rameter ([13]) is given by z max = 2ZZ^e 2 R/[Uvb m i n ('yo + 
71)] with the target charge number Z — 82 and the effec- 
tive charge number — 4/11. The numerical values for 
•Zmax are given again in table [JJ They are smaller than one 
and decrease with the projectile energy. The range where 
the reduction factor is actually needed is well approxi- 
mated by the NLO correction, which corresponds to the 
approximation of r(z) in figure[TJby an inverted parabola. 
In the present approach we are more certain about the 
nuclear structure input that determines the opening pa- 
rameter of this parabola. 

Similar as in [5] we define a reduction factor R(v) = 
tj(°°) / u LO of the total excitation cross section er(°°) caused 
by higher-order effects. This factor R(v) can be used to 
take higher order effects into account in an analysis of the 
experimental data. 

Since the parameter z is always smaller than one in 
the experiments, the full cross section in sudden approxi- 
mation ct' 00 ' = <jlo + (?nlo + • ■ • is well approximated by 
the NLO correction and the reduction factor becomes 



R(v) « 1 + 



1 



(26) 



The numerical values as given in table [TJ show that the 
reduction of the first-order cross section is in the order 
of 6 to 11 % for the various experiments. In figure [2] the 
dependence of the reduction on the projectile velocity v 
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Fig. 2. Reduction factor R(v) = 1 + onlo /<tlo as a function 
of the velocity v for values of (dashed line), 0.604 (solid line) 
and 1 (dotted line) of the ratio 71 /70 = \J ' E\/Eq. 



is depicted for the actual value of the ratio 71 /70 and the 
two limiting cases as in figure [1] The curves are shown 
only for velocities where the parameter z does not exceed 
the value one, the radius of convergence for the expansion 
of the amplitude (jT5]> . Note that an increase of 71 from 
zero (dashed line) to one (dotted line) leads to a reduction 
of the parameter z and correspondingly to a reduction of 
the higher-order corrections in the cross section for the 
same velocity v. 



4 Conclusion 

We have presented a realistic model for higher-order ef- 
fects in the Coulomb excitation of neutron halo nuclei. 
With quite simple methods a reliable value for the re- 
duction of the cross section could be given. In addition, 
this should be very useful as a benchmark, for valuable 
and necessary tests of more sophisticated approaches like 
CDCC or time-dependent calculations. These kinds of cal- 
culations require a lot of computation, yet they should 
reproduce the present results in certain well defined lim- 
its. Using an effective-range point of view we updated our 
previous work [5] . Other angular momentum combinations 
can be treated in a similar way. However, for higher an- 
gular momenta the halo nature is less pronounced and 
the importance of the core size parameter R increases. 
Higher-order terms in 7 will become more important and 
the method less useful. It is quite remarkable that our 
results for 11 Be Coulomb excitation depend only on the 
binding energies (they determine the parameters qo and 
<?i) and a parameter R which characterizes the core size. 
This stresses the fact that halo nuclei depend only on a 
few low-energy constants, and not on details of the shape 
of the potential. This is in accord with low-energy scat- 
tering, which is also characterized by a few low-energy 
parameters in the effective-range expansion. 

In view of the discovery of the neutron-rich isotopes 
40 Mg and 42 Al [2TJ1I2T] one may expect that more neutron- 



6 



S. Typel, G. Baur: Scaling laws and higher-order effects in Coulomb excitation of neutron halo nuclei 



rich nuclei with low-lying single particle states exist to 
which the present approach can be applied. It is very likely 
that the new radioactive beam facilities around the world 
like FAIR, SPIRAL2 and RIBF will provide interesting 
new examples in the higher mass region. 
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Appendix 

For the analytical calculation of the excitation amplitude 
in sudden approximation the argument of the exponential 
function in (jHJ) can be written as qcoui ■ T = qcouir cos(#) 
in spherical coordinates because the vector qcoui points in 
z-direction. Using 

yi f f s mf (r)]' y^ mi (r) = (-iy^5 mtmf C -^f> (27) 



with a spherical Bessel function j\ (x) — — j^jo (x) so that 



We notice 



lim I 2 (p,q) = 0. 

q^O 



-i-h(p,q) = --h(p,q) 
dp q 



(33) 
(34) 



and find with the help of a partial integration 



h(p,q) = C--[dp , I 1 (p',q) (35) 
1 Jo 

= C — 2- + arctan ( - 



q q 

H arctan - 



q \q 

where a constant C appears. Since 



arctan(x) + arctan I — ] = — 
x) 2 



(36) 



for = j f = \mi\ = \m f \ = \, k = 0, l f = 1 and s = 1/2, for a; > we can write 
one finds 

a fi 

= (-1) 



(28) 



ji+m 



CqCi 

47T 



OO />2"7T />1 

dr I d(j) I dcos(9) cos(0) 



o 



hM = C-* q+ (l + ^\ arctan (J j - | . (37) 



With 



x exp[-(i?i + q )r] 1 H exp[-iq Cou ir cos(6»)] 

q x r ; 



\imI 2 (p,q) = C--=Q 



(38) 



CoC\ 
2i 



1 1 
h + -h 

qcon\ qi 



with two integrals that can be calculated separately. The 
first integral is given by 



the constant can be determined to be C = tt/2 and the 
final result is 



I, [p. <j) j 1 + arctan f |J ~ | ■ (39) 



p OO pi 

h(p,q) = iqj dr J dz z exp[-(p + iqz)r] (29) References 



iq I dz 

i p + iqz 

iq p — iq 

with p = go + qi and q — gcoui ■ With the relation 
1 + iz 



In- 



1 



iz 



2i arctan(z) 



one obtains 



hip, q) 



P ,11 
— arctan 

q 



The second integral is given by 



h (p, q) = i J dr J dz z 



'' exp[— (p + iqz)r] 



dq Jo 



1 exp[— (p + iqz)r] 



= 2 



dr jiiqr) 
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